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Abstract. We show how a difference in the cor- 
relation length longitudinally and transversely, 
with respect to the jet axis in e + e~ annihilation, 
arises naturally in a model for Bose-Einstein cor- 
relations based on the Lund string model. The dif- 
ference is more apparent in genuine three-particle 
correlations and they are therefore a good probe 
for the longitudinal stretching of the string field. 

1 Introduction 

Bosons obey Bose-Einstein statistics, which 
compared to an uncorrelated production leads 
to an enhancement of the production of iden- 
tical bosons at small momentum separations. 
The enhancement is called the Bose-Einstein 
effect and it is very often parametrised in the 
form 

R 2 (q) = l + Xcxp(-Q 2 R 2 ) (1) 

where Q is the relative four-momenta of a pair 
of bosons, Q 2 = —q 2 = — {P2—P1) 2 , arid R and 
A are two phenomenological parameters. The 
parameter R is often refered to as the radius 
of the boson emitting source. It is clear that 
the correlation function R 2 reflects the space- 
time region in which the particle production 
occurs but both the paramctrisation in equa- 
tion 1 and the interpretation of R are often 
given without very convincing arguments. 

In reference [1] (an extension of reference [2] 
to multi-boson final states) a model for Bose- 
Einstein correlations based on a quantum- 
mechanical interpretation of the Lund string 
fragmentation model [3] is presented. In this 
work we will investigate some features of the 



model to see how the correlation lengths in 
the model arises and this will be used to show 
what the parameter R is sensitive to. We will 
in particular show that the model predicts, 
due to the properties of string fragmentation, 
a difference between the correlation length 
along the string and transverse to it. In prac- 
tice this means that if we introduce the longi- 
tudinal and transverse components of the vec- 
tor q (defined with respect to the thrust axis in 
e + e _ annihilation) we obtain a noticablc dif- 
ference in the correlation distributions. This 
becomes even more apparent when we go to 
three-particle correlations because in this case 
one is even more sensitive to the longitudinal 
stretching of the string field. 

2 Correlation lengths 

The starting point of our Bose-Einstein model 
[1, 2] is an interpretation of the (non- 
normalised) Lund string area fragmentation 
probability for an n-particlc state (see fig- 
ure 1) 

dP(p 1 ,p 2 ,...,Pn) = (2) 

n 

Y\Nd Pj 6 (p 2 -m 2 )5{Y,P 3 - Ptot )exp{- bA) 
i 

in accordance with a quantum mechanical 
transition probability containing the final 
state phase space multiplied with the square 
of a matrix element A4. In reference [2] and 
in more detail in reference [1] a possible ma- 
trix clement is suggested in agreement with 
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(Schwingcr) tunneling and the (Wilson) loop 
operators necessary to ensure gauge invari- 
ance. The matrix clement is 



M = expire - b/2)A 



(3) 



where the area A is interpreted in coordinate 
space, k ~ 1 GeV/fm is the string constant 
and b ~ 0.3 GeV/fm is the decay constant. 

The transverse momentum properties are 
in the Lund model taken into account by 
means of a Gaussian tunneling process. The 
produced (qq)-pair in each vertex will in this 
way obtain ±k^ and the hadron stemming 
from the combination of a q from one ver- 
tex and a q from the adjacent vertex obtains 

P-L = k ^j + l - k -Lj- 

In case there are two or more identical 
bosons the matrix element should be sym- 
metrised and in general we obtain the sym- 
metrised production amplitude 



(4) 



where the sum goes over all possible permuta- 
tions of identical particles. Taking the square 
we get 



\Mf = J2\Mv\ 2 1+E 



■p'jiV 



2Rc(M v M* v ,) 
\Mv\ 2 + \M v ,\ 2 

(5) 

The MC program JETSET [4] will provide the 
outer sum in equation 5 by the generation of 
many events but it is evident that the model 
predicts a quantum mechanical interference 
weight, w-p, for each given final state charac- 
terised by the permutation V: 



W-p 



1 



E 

pi^p 



2Re(M v M v ,) 
\Mp\ 2 + \M V >\ 2 



(6) 



In the Lund model wc note in particular for 
the case exhibited in figure 1, with two iden- 
tical bosons denoted 1 and 2 having a state / 
in between, that the decay area is different if 



the two identical particles are exchanged. It is 
evident that the interference between the two 
permutation matrices will contain the area dif- 
ference, AA, and the resulting general weight 
formula will be 

AA 



Wp 



cos ■ 



-P'^-P cosh 



bAA 



+ 



2H 



(7) 

where A stands for the difference between the 
configurations described by the permutations 
V and V' and the sum is taken over all the ver- 
tices. The calculation of the weight function 




Fig. 1. The. decay of a Lund Model string span- 
ning the space-time area A. The particles 1 and 
2 are identical bosons and the particle(s) produced 
in between them is denoted by I. The two ways to 
produce the state are shown. 



for n identical bosons contains n! — 1 terms 
and it is therefore from a computational point 
of view of exponential-type. We have in refer- 
ence [1] introduced approximate methods so 
that it becomes of power-type instead and we 
refer for details to this work. 

We have seen that the transverse and longi- 
tudinal components of the particles momenta 
stem from different generation mechanisms. 
This is clearly manifested in the weight in 
equation 7 where they give different contri- 
butions. In the following we will therefore in 
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some detail analyse the impact of this differ- 
ence on the transverse and longitudinal corre- 
lation lengths, as implemented in the model. 

In order to understand the properties of the 
weight in equation 7 we again consider the 
simple case in figure 1. The area difference of 
the two configurations depends upon the en- 
ergy momentum vectors p\ , p 2 and pi and can 
in a dimensionless and useful way be written 

AA 

^ = 5 P 5x L (8) 

where 5p = P2 — Pi and Sxl — (St;0,0,8z) 
is a reasonable estimate of the space-time dif- 
ference, along the surface area, between the 
production points of the two identical bosons. 

To preserve the transverse momenta of the 
particles in the state (1, /, 2) it is necessary to 
change the generated at the two internal 
vertices around the state I during the permu- 
tation, i.e. to change the Gaussian weights. 
Also in this case we may write a formula simi- 
lar to equation 8 for the transverse momentum 
change: 



2k 



— Sp±5x± 



(9) 



where Sp ± is the difference p±2 — Pxi and 
Sx± = (k±b — (— kj_ ))/«. The two neighbour 
vertices to the state (1, I, 2) ((2, /, 1)) are de- 
noted by a and b and kj_b+k^ a corresponds to 
the state's transverse momentum exchange to 
the outside. Therefore 5x± constitutes a pos- 
sible estimate of the transverse distance be- 
tween the production points of the pair. 

For the general case when the permutation 
V' is more than a two-particle exchange there 
are formulas similar to equations 8 and 9. 

It is evident from the considerations leading 
to equations 8 and 9 that only particles with 
a finite longitudinal distance and small rela- 
tive energy momenta will give significant con- 
tributions to the weights. We also note that 
we are in this way describing longitudinal cor- 
relation lengths along the color fields, inside 



which a given flavor combination is compen- 
sated. The corresponding transverse correla- 
tion length describes the tunneling (and in 
this model it provides a damping chaoticity). 

The weight distribution we obtain is dis- 
cussed in reference [1]. It is strongly centered 
around unity and we find negligible changes in 
the JETSET default observables (besides the 
correlation functions) by this extension of the 
Lund model. 

3 Results 

We have analysed two- and three-particle cor- 
relations in the Longitudinal Centre-of-Mass 
System (LCMS). For each pair (triplet) of par- 
ticles the LCMS is the system in which the 
sum of the particles momentum components 
along the jet axis is zero. In the pair analysis 
we have used the kincmatical variables 



q± = \J (Px2 - Pxl) 2 + (Py2 - Pyl) 2 (10) 
<7L = 1^22 ~Pzl\ 

and in the triplet analysis we have used 



9x 



9±12 + 9±13 + 1±23 (11) 



where the jet axis is along the z-axis. 

We have taken the ratio of the two-particle 
probability density of pions, p 2 , with and 
without BE weights applied as the two- 
particle correlation function, R2 



R2{Pl,P2) 



P2w(j>\,P2) 
P2{Pl,P2) 



(12) 



and the resulting distribution is shown in fig- 
ure 2. It is clearly seen that it is not symmetric 
in qz, and q± and in particular that the corre- 
lation length , as measured by the inverse of 
the width of the correlation function, is longer 
in the longitudinal than in the transverse di- 
rection. This difference remains for reasonable 
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changes of the width in the transverse mo- 
mentum generation. Using all the final pion 
pairs in the analysis results in in a small de- 
crease in the transverse correlation length and 
of course a large decrease in the height for 
q L ~ q± ~ 0, while the longitudinal correla- 
tion length is rather unaffected. 



i? 3 = i?3 - i?3 + f (15) 

This way of getting the genuine correlations 
is not possible in an experimental situation, 
where one has to find other ways to get a 
R 3 reference sample. We have suggested one 
possible option in reference [1]. The distribu- 
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Fig. 2. The two-particle correlation function 
R,2{qL,q±) for charged pions. The sample con- 
sistes of particles which are either initially pro- 
duced or stemming from short-lived resonances. 



The total three-particle correlation function 
is in analogy with equation 12 



# 3 (Pl>P2,P3) 



P3m(Pl,P2,P3) 
P3(Pl,P2,P3) 



(13) 



To get the genuine three-particle correlation 
function, R 3 , the consequences of having two- 
particle correlations in the model have to be 
subtracted from i? 3 . To this aim we have cal- 
culated the weight taking into account only 
configurations where pairs are exchanged, w . 
In this way the three-particle correlations 
which only are a consequence of lower order 
correlations can be defined as 



R 3 (Pl,P2,P3) 



P3(Pl,P2,P3) 



(14) 



The genuine three-particle correlation func- 
tion, i? 3 , is then given by 
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Fig. 3. The genuine three-particle correlation 
function Rs(qL,q±) for all final state charged pi- 
ons. 



tion i? 3 is shown in figure 3. The effect of the 
higher order terms is to pull the triplets closer 
in longitudinal direction while the transverse 
direction is rather unaffected. This suggests 
that the higher order terms are more sensi- 
tive to the longitudinal stretching of the string 
field. 
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